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ABSTRACT 

Rigcl (/3 Ori, B8 la) is a nearby blue supcrgiant displaying a Cyg type variability, and is one of 
the nearest type-II supernova progenitors. As such it is an excellent test bed to study the internal 
structure of pre core-collapse stars. In this study, for the first time, we present 28 days of high precision 
MOST photometry and over 6 years of spectroscopic monitoring. We report nineteen significant 
pulsation modes of SNR>4.6 from radial velocities, with variability time scales ranging from 1.21 to 
74.7 days, which are associated with high order low degree gravity modes. While the radial velocity 
variations show a degree of correlation with the flux changes, there is no clear interplay between the 
equivalent widths of different metallic and Ha lines. 
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1. INTRODUCTION 

Blue supcrgiants (BSGs) such as Rigcl are post-main 
sequence (MS) massive stars (M > 15M©) evolving 
quickly across the Hcrtzsprung-Russcl diagram. They 
are amongst the most intens ively studied objects in c on- 
temporary astrophysics (see lMaeder fc MevnetH20iTl for 
a recent review). BSG stars are intriguing objects be- 
cause they end their lives as Type II supernovae (SN-II) 
which serve as catalysts for star formation and primary 
contributors to the heating and chemical enrichment 
of their host galaxies. These stars also offer a poten- 
tially valuable addition to the distance l adder by way of 
the B SG wind-momentum-luminosity dKudritzki et al.l 
I1999D and flux-weighte d gravity-luminosity relationships 
(jKudritzki et al.l 120031) which give very consistent dis- 
tance determinations to the Local Group of galaxies 
dKudritzki et all 120081: III et alJl200lh . 

Fascinating featur es of BSGs in clude their microvari- 
ability in flux (e.g. iSterkenl Il977l ). complex changes in 
radial velocity (RV), as well as variations in the equiv- 
alent widths (EW) and line profiles of Balmer (espe- 
cially Ha) and metallic lines (jKaufer et all I1996L ll997t 
iRichardson et all 120111 ) . Kaufer and coworkers moni- 
tored the optical spectra of six BA supcrgiants (includ- 
ing Rigcl and a Cyg). They analyzed the radial ve- 
locity measures of these stars and associated their ob- 
served variability with non-radial pulsations since the 
traveling features in the dynamical spectra could not 
be reconciled with th e rotational periods of the stars. 
IWaelkens et al.1 (|1998| ) discovered 32 pulsating BSGs in 
the Hipparcos database which belong to the a Cyg class 
of variable stars, with periods ranging from 1.5 to 24 
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days. They were later definitely identified as gravity 
mode oscillations by a Non-LTE spectroscopic analysis of 
iLefever et al.l (|2007|) who compared the current position 
of their stars with the instability strip of g-mode domi- 
nated pulsators. Beyond the Milky Way, iBresolin et al.l 
(|2?)0l detected a handful of variable BSGs in NGC 300, 
and found the perio ds of two of these star s to be 72.5 and 
96 days. Recently, lAerts et "all (|2010b( ) correlated the 
sudden amplitude decline in the spectroscopically pecu- 
liar CoRoT target HD 50064 (Bl-6 la) to radial strange 
mode variability with a 37 day period; they inferred a 
mass of ~ 45A/ Q for this luminous mass-losing BSG. De- 
spite several decades of observational efforts on ground 
based photometry and spectroscopy of a Cyg type pul- 
sating stars, it is not certain what portion of their vari- 
ability is periodic, nor how far they deviate from strict 
periodicity. 

The nature of the aforementioned variability in high 
mass post-MS variables is still poorly understood. The 
photometric variability of HD 163899 (B2 lb /ll) found 
by MOST observations reve aled a total of 48 pressure 
(p) and gravity (g) modes. iSaio et all (|2006[ ) showed 
that these modes can be simultaneously excited by the 
K-mechanism and reach the surface if they arrive at the 
base of the intermed iate conve c tion z one (ICZ) with an 
appropriate phase. iGautschvl (|2009f ) searched for the 
origin of long-period variabilities in the prototype a 
Cyg. Interestingly his Figure 5 shows a gap where no 
instabilities are predicted for evolutionary tracks with 
3.95 < logT e ff < 4.15. Rigel lies in this gap. In con- 
trast to this result, the study of lSaiol (|2011[) predicts that 
Rigel should be unstable against non-radial convective 
g - modes. This contradiction could arise from different 
physical ingredients (such as rotational and overshoot 
mixing) and various numerical techniques in pulsation 
codes. This is not surprising since evolved massive stars 
like Rigel arc demanding to model. 

The radial strange mode is proposed as another mech- 
anism to induce instability and interplay with mass 
loss efficiencies in these massive stars with log(L/M) > 
4 log( W M^) (see e.g ~ 
ISaio! 120 lit lAerts et al.1 
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2010bL and references therein) 



This requirement is also marginally fulfilled by Rigel. 
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IGodart et al.l (|2009( ) investigated the destructive impact 
of core overshooting and mass loss during the MS life- 
time on the extent of the ICZ, and showed that models 
with wider ICZs are more likely to destabilize stellar os- 
cillations. These studies show that asteroseismology of 
slowly pulsating blue supergiants can reveal a wealth of 
information about the internal structure of SN-II progen- 
itors. This can be complementary to the understanding 
of the physical properties o f pulsars which ori ginate from 
the cores of massive stars (jHeger et al.l 120051 ). 

This is the first paper in a series of investigations aim- 
ing at probing the details of the internal structure of 
BSGs through asteroseismic study, and investigating the 
origin of their short- and long-period variability. Rigcl 
was selected as an ideal test subject because of its appar- 
ent brightness, proximity and current evolutionary sta- 
tus. The latest measurements of the physical parame- 
ters of Rigel collected from the literature are summarized 
in Section [5J The space-based MOST photometry and 
ground-based spectroscopy are presented in Section [3l 
The results of the multimode pulsation frequencies are 
finally discussed in Section [4] The interpretation of the 
pulsation frequencies will appear in Moravveji et al. (in 
preparation) . 

2. ABOUT RIGEL 

Rigel (fi Ori; HD 34085; B8 la; V-0.12 mag) is the 6 th 
brightest star in the night sky and the most luminous star 
in the solar neighborhood. It is a member of a multiple 
star system where its companion, Rigel B. is a spectro- 
scopic binary about 9.5"distant (|Sanfordlll942[ ). From 
the revised Hipparc os parallax the dis tance to Rigel is 
dmp = 264 ± 24 pc (|van Leeuwenll2007fl which is s maller 
than 360 ± 40 pc adopted bv iHofHeit fc Jaschekl (|1982l ) 
who assumed that Rigcl originates from the r Ori Rl 
complex. 

Fundamental parameters of Rigcl which impose valu- 
able constraints on the equilibrium model of the star 
are already measured. They are effective temperature 
T cff = 12 100 ±150 K, surface gravity log g = 1.75 ±0.10, 
luminosity log(L/L0) = 5.08+q ioi near-solar metalic- 
ity [M/H]=-0.06 ± 0.10, surface He a bundance Y s = 
0.32 ± 0.04 vsini » 25 ± 3 km s" 1 (jPrzvbilla et al.l 
[200l [20101: iSimon-Diaz eFail 1201(1 . The most up-to- 
date limb darkened angular diameter for Rigcl comes 
from CHARA/ FLOUR K-band interfe rometry 6> LD = 
2.75 ±0.01 mas (jAufdenberg et alJ[2008l ) which in combi- 
nation with dmp yields R = 78.9±7.4i? Q . iPrzvbilla et al.1 
propose a mass of 23Mq from Geneva evolution- 
ary tracks which include the effects of rotation. 

The absolute bolometric magnitude is ca lculated 
(adop ting the observed B — V = —0.03 from iNicoletl 
(11978[l. and a bolomet ric correction BC = —0.78 
from iBessell et al.1 (fl998h with E(B - V) = 0.05 from 
IPrzvbilla et al.l ([2006^ to yield M bo i = -7.92 ±0.28. By 
applying the flux-weig hted gravity-luminosity relation of 
iKudritzki et~all (|2003l) 

M bol = 3.71 log(. 9 /T c 4 ffi 4 ) - 13.49 

(T eff , 4 = T cff /10 000) which is optimized for BSGs, we 
arrive at a value of Mboi = —8.17 ± 0.45 - in good agree- 
ment with our calculated value. 
Moreover, Rigel's line profile variability in Ha 



TABLE 1 

Updates to Physical Parameters of Rigel 



Parameter Value Reference 



"Hip lP c l 
T cS [K] 


Z04 ± Z4 


van Lccuwcn (2007) 


12 100 ± 150 


Przvbilla et al. (2010) 


]oe:(L/Lr.) 
iU &V- L '/ -^0 J 


J ' uo — 0.10 


This work 


Y s 


0.32 ± 0.04 


Przvbilla et al. (2010') 


v sin i [km s ] 


25 ±3 


Przvbilla et al. (2010) 
Simon-Diaz et al. (2010) 


Inclination 


60° < i < 90° 


This work 


6» LD [mas] 


2.75 ±0.01 


Aufdenbere et al. (2008) 


R/Rq 


78.9 ± 7.4 


This work 


M bol I ma g] 


-7.92 ±0.28 


This work 


M [M Q yr" 1 ] 


l-2xl0~ 7 


Chesneau et al. (2010) 


B d [G] 


< 25 


Shultz et al. (20lTl 


-^shortest [d] 


1.2191 ±0.0001 


This work 


-^longest M 


74.74 ± 0.28 


This work 



([Kaufer et al.l I1996L iMorrison et~aT1 I2008Q and other 
metal-line equivalen t widths due t o non-radial pulsations 
(|Kaufer et al.lll997h and mass loss ()Chesneau et al.|[2010D 
is already published. 

Rigcl shows variations in the Ha spectr al feature with 
some outburst events being recorded bv llsraelian et al.l 
(|1997|) and Chesneau (VLTI campaign, i n preparation). 
The m orphology of this line is studied bv [Morrison et al.l 
(|2008l) . Moreover, iChesneau et al.l (|2010D detected an 
extended rotating Ha r egion from optical in terferometry 
with VEGA/CHARA (jMourard et al.l [2009]) and specu- 
late that Rigcl is observed with its rotation axis along 
the North-South direction at significantly high inclina- 
tion angle. They estimated the current mass loss rate 
is M = 1 — 2 x 10~ 7 M© yr -1 . Assuming a high incli- 
nation angle, the rotation velocity is well below 20% of 
the estimated cri tical breakup r otation rate of ~185 km 
s _1 (Eq. 2.19 in [Macdcr 2009|), so spherical symmetry 
can be safely assumed for its geometry. The current best 
values of the physical parameters of Rigcl are presented 
in Table [U Additiona l parameters for Rigel are given by 
IPrzvbilla et al.l (|2006D . These are invaluable input for 
future seismic modeling and analysis. 

3. OBSERVATIONS 

3.1. MOST High Precision Photometry 

Rigel was observed continously with the MOST satel- 
lite (Microvariabi lity and Oscillations of Stars, 
iWalker et all l200l for 27.7 days from 15 Novem- 
ber to 13 December 2009. The dataset consists of 
30640 observations after correcting for the Southern 
Atlantic Anomaly. Figure [T]a shows the light curve 
after removing an offset of ~ 0.12 mag. The abscissa 
is in HJD-2 450 000.0. Because of the brightness of 
the target, the precision of the observations is of the 
order of «0.05 to 0.10 mmag. Starting with a sinusoidal 
variability pattern, oscillations die out and Rigel seems 
quiescent for several days, approximately from HJD 
5164 to 5168 (enlarged in Figure [TJb), followed by a 
gradual rise, and a steep decline in light flux; this 
resembles a possible beating pattern and we expect 
close, low-frequency modes to appear in the harmonic 
analysis (Section^]). There is a 20-hour gap around HJD 
5175, arising from an interruption in communication 
between the satellite and the ground station. This is 
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Fig. 1. — Compilation of observations from the Rigel campaign. Panel (a) shows the 27.7-day MOST space photometry (Section I3.lt . 
Panel (b) is a 5-day zoomed-in portion of the data to emphasize the amplitud e and time scale of the shortest variations. Panel (c) shows 
the RV data which demonstrates a correlation with the flux changes ( Section 13.21 1. Panel (d) is the equivalent width variation for Si II 
A6347.1 A during the same monitoring epoch (Section l3,4t : the (red) squares are the 1-day binned averages (sec Figure 



at the time when the star is decreasing in brightness 
almost monotonically. 

3.2. Spectroscopy 

Rigel's optical spectrum has been monitored by the 2- 
m Tennessee State University Automatic Spectroscopic 
Telescope at Fairborn Obse rvatory, Arizona (TSU:AST, 
lEaton fc Williamson! I2007I ). A total of 2 328 high- 
resolution (R~30 000 and 20 000) moderate signal-to- 
noise (SNR-50 to 150) echelle spectra (4 900 - 7100 A) 
were obtained. The spectra were secured over 334 nights 
starting from 11 December 2003 to 14 February 2010. 
The density of time-sampling was increased during the 
months centered on the MOST observations. Simulta- 
neous with MOST, 442 spectra were collected over 20 
nights. 

The spectra were reduced at TSU by an automatic 
pipeline method without removing the telluric lines. The 
radial velocity variations were derived by fitting a Gaus- 
sian to 29 metallic lines; Figure [TJc shows the RV vari- 
ations contemporaneous with MOST observations. An 
offset of +0.3 km s _1 w as applied to the AST instr umen- 
tal RV measurements (lEaton fc Williamson] I2007f ). The 
standard deviation in the mean for each RV measurement 
was evaluated as the standard deviation of the 29 mea- 
sured absorption lines in each spectrum. Averaged over 
all spectra, this value is 0.26 km s _1 . In the years prior 
to the MOST observations, Rigel was observed sporadi- 
cally with AST. Also, the AST dataset is interrupted by 



cloudy weather that occurred during the MOST observ- 
ing run. Thus the sampling is not even, and the analysis 
of these data will suffer from daily and annual aliasing. 

Extensive and highly variable velocity changes are 
clearly visible in Figure [TJ:. During the plotted interval, 
the RVs vary from 11 to 25 km s . There is evidence for 
a correlation between the RV and MOST brightness vari- 
ations. During this interval, and other seasons, the RV 
variations appear to exhibit complex periodicity typical 
of a multi-periodic pulsating star. 

In addition to the AST spectroscopy, Rigel was ob- 
served with the VLTI during the time of the MOST ob- 
servations. A preliminary analysis of the spectro- 
interferometry by Chesneau et al. (in preparation) shows 
that during this interval Rigel was relatively quiet, with 
the Ha line appearing in absorption with sometimes a 
weak emission component appearing at times, and no 
outbursts being recorded. No major high velocity ab- 
sorption (HVA) event is identified during this interval. 
This is confirmed by seven raw Ritter echelle spectra (R 
= 26,000, SNR 100 to 200) taken during the same inter- 
val at Ritter Observatory (Morrison; private communi- 
cations). As a result, there is no evidence of a possible 
propagating atmospheric shock or chaotic activity during 
MOST photometry. 

3.3. Spectropolarimetry 

In conjunction with MOST photometry, 78 observa- 
tions were collected with the high-resolution (R-~65 000) 
spectropolarimeter ESPaDOnS, installed at CFHT, and 
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Fig. 2. — (a) A representative 50 day radial velocity curve of 
Rigel. 1-ct error in RV is ovcrplotted. (b) A residual plot with 19 
harmonics removed as in Tablc[2]with rms deviation 0.59 km s _1 . 
(c) The 1-day binned EW measures of two metallic lines, namely 
Si II A 6347 and Si II A 6371 A which exhibit significant temporal 
variations, (d) A similar plot for the He I A 6678 line, and (e) for 
C II A 6578 and Ha lines. 

its clone Narval at Telescope Bernard Lyot (TBL, 
Pic du Midi). Analysis with the multi-line cross- 
correlation techniq ue Least Squares Deconvolution (LSD; 
IDonati et al.lll997D reveals no evidence of a magnetic field 



( Shultz et al.l 1201 If ) . with error bars on the longitudinal 
field B[ of order 15 G. Matching synthetic disk-integrated 
Stokes V profiles to the observed Stokes V profiles has 
constrained the dipolar magnetic field to be Bd < 25 G 
for for high inclinations of the rotational axis and low 
obliquities of the magnetic axis, although Bd ~ 50 G is 



possible at low-intermediate inclinations, in which case 
the rotational period is shorter and fewer observations 
can be binned together (Shultz et al., in preparation). 
If, as discussed above, Rigel is viewed nearly equator on, 
then the first of the conditions necessary for a higher 
upper limit applies. However, these upper limits cannot 
rule out the possibility of significant magnetic effects. 
Thus, in the absence of any positive evidence for a mag- 
netic field, or any pressing theoretical reason to suspect 
its existence, there is no reason at this time to complicate 
the pulsational analysis with its inclusion. 

3.4. Equivalent Widths 

iKaufer et al.1 (|1997|) showed that different spectro- 
scopic lines in BA supergiants emerge from different op- 
tical depths t\ , in a sense that lines from shallower pho- 
tospheric depths have larger equivalent widths (hereafter 
EW designated by W\). They classified the lines in the 
optical band accordingly to weak (50 < W\ < 200 mi), 
medium (200 < W\ < 500 ml) and strong (W x > 500 
mA). We chose C II A6583 as a weak line, Ha, C II 
A6578 and Si II A6371 as medium lines, and He I A6678 
and Si II A6347 as strong lines. 

EWs of the above lines are calculated according to 



W x = /(l — F\) dX. 



(1) 



where F\ is the flux at each wavelength interval dA, 
renormalized to the continuum on either side of the spec- 
tral line. To suppress the contribution from cosmic rays 
the spectra were filtered using a low-bandpass median 
filter; contamination from telluric lines was removed by 
identifying those regions so contaminated and using a 
higher-bandpass median filter. The corresponding errors 
for each EW measurement are evaluated with the per- 
pixel standard deviation of the flux from the mean value 
of a nearby section of the continuum and propagated 
through the measurement. Uncertainty in the location 
of the continuum is accounted for in the calculation of 
the per-pixel flux error from the SNR. As an example, 
Figure QJ1 shows the variations in W\ for the Si II A6347 
A line during the MOST photometry. There is a weak 
correlation between the RV and EW variations during 
this season but it is not always present. 

Comparison of ESPaDOnS/Narval spectra with con- 
temporaneous AST spectra reveal that, while the instru- 
ments generally yield consistent measurements of W\, 
the comparatively lower-resolution AST spectra show a 
small but persistent bias in which W\ systematically in- 
creases with the uncertainty a(W\). However, W\ is 
accurate to within ~5% which is sufficient to reveal at 
times a weak correlation between RV and the W\ of, for 
instance, the Si II 6347 A line, one of the stronger pho- 
tosphcric lines in the AST spectral window (see Fig. [T}1) . 
Figure[2]is a 50-day representative interval, occurring one 
year prior to the MOST photometry, that shows the RV 
changes as a function of time (panel a). For the same 
time interval, panels (c) to (e) plot the equivalent widths 
of five different lines. While Ha (having high sensitivity 
to wind) demonstrates the largest variation amplitude, 
the rest of the lines, though they have different absolute 
amplitudes, change moderately - on the order of 12 to 
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19%. The two Si II lines follow a similar trend and have 
similar amplitudes. 

4. RIGEL PERIODICITIES 
4.1. Variability Pattern of Rig el 

iKaufer et al.l (|1996l I1997D derived several harmonics 
to fit to the season-by-season photometry, radial veloc- 
ity and EW variations of Rigcl and five other late BA 
supergiants. As they demonstrated, the Ha line equiv- 
alent widths show systematic variations. Figure 7 in 
IKaufer et al.l (119971 ) shows that a Cyg, the prototype of 
this class of pulsating stars, and the other BSGs exhibit 
variability with periods exceeding a week. For the spe- 
cific case of Rigel, the per iods they find are in the r ange of 
4 to more than 50 days. iRichardson et all (|2011h found 
evidence for season-to-season changes in the oscillation 
frequencies of a Cyg (Deneb, A2 la) from five years of 
spectroscopic and photometric monitoring. 

The need for multiple modes to fit to our RV and 
EW data set is clearly shown in Figure [2] Similar be- 
havior is observed in other seasons. Because different 
spectral lines are formed at various optical depths of the 
star's atmosphere, the temporal changes in W\ for strong 
lines (such as Si II A6347) are significantly different from 
those of weak lines (such as C II A6578). The equiva- 
lent width, however, is sensitive to local changes in the 
effective temperature ST c ff/T e g and gravity 5g/g. In a 
series of studies in (3 Cep and SPB stars (which are dif- 
ferent from BA supergiants, b ut still have similarit ies in 
the nature of th e ir pul sation) , iDupret et al.l (|2002f) and 
De Riddc r~et al.1 (|2002t ) showed that a sinusoidal behav- 
ior in EW and RV time series with a common frequency 
is observed in only few cases. Given that Rigel (similar 
to SPBs) pulsates in g-modes which are transversal mo- 
tions at t he stellar surface, EWs are strongly affected by 
ST cS /T eS (|De Ridder et alj|2002ft . Based on the observed 
behavior of Rigel's equivalent widths (panels c to e in Fig- 
ure^, we believe that a complex non-r adial velocity field 
exists across different o ptic al depths (lAerts et al.1 120091 : 
iSimon-Di'az et al.ll2010l ). As IKaufer et all (|1997| ) showed 
(in their Figure 4), radial velocity variations are negligi- 
bly affected by this depth dependence. Furthermore, the 
RV variations show a weak (not to say absent) degree of 
correlation with the EW of any of the lines. 

4.2. Frequency Analysis Based on RV data 

Compared to the MOST photometry time series 
(Rayleigh limit of 0.036 d -1 ), the RV data have a much 
longer time span (Rayleigh limit of 4.43 x 10~ 4 d _1 ). 
As a result, we base our search for intrinsic signals pri- 
marily on this data set usi ng two widely -used state-of- 
the-art programs - Si gSpec (Recgen 2007) and Period04 
(|Lenz fe Br cger 2005|). However, the RV data suffer from 
strong daily (fd = 1 d _1 ) and annual (f a = 0.002769 
d _1 ) aliasing (inner panel in Figure [3]). As shown in 
Figure [3l there is a repeating pattern every 5 d _1 in 
the discrete Fourier transform (DFT) spectra that arises 
from our sampling rate. Thus we selected the upper fre- 
quency scan range at / nyq =5 d , and our analysis did 
not identify frequencies higher than Id -1 . 

SigSpec computes the spectral significance level (sig) 
for the DFT amplitude spectrum based on the analyti- 
cal solution for the Probability Density Function (PDF) 




Frequency [cd"'] 

Fig. 3. — Discrete Fourier Transform of the RV time series shows 
a repeating pattern at multiples of /nyq = 5 d . The inner panel 
shows the spectral window and the one day aliasing in addition to 
the annual aliasing. 




Frequency [cf 1 ] 

Fig. 4. — Significance spectrum of RV data with a normalized 
weight. Only the low-frequency domain between [0, 1] d _1 is 
shown. The inner panel shows the remaining power below the 
threshold sig=15. After prewhitening of 19 significant harmonics, 
the residual is 0.94 km s — . 

of an amplitude level of any peak during the iterative 
prewhitening process. By default the program scans for 
peaks with sig> 5 which is approximately equivalent to 
S/N= 3.8 ( Rccgcn 120071 ) . However, we conservatively 
chose just to search for highly significant modes to avoid 
a forest of low-frequency peaks. To accomplish this, the 
prewhitening process is stopped if the significance (sig) 
of each peak or the cumulative significance (csig) of the 
wh ole solution is be l ow 15 . A similar threshold was used 
by iChapellier et al.l (|2011l ) in the frequency analysis of 
one of the CoRoT primary targets. Since each RV mea- 
surement i has an error o~i, we associate a normalized 
weight Wi = a^ 2 / er~ 2 to each data point. The result- 
ing significance spectrum is shown in Figure 2) This pro- 
cedure leads to a prewhitening of 19 significant modes. 
The corresponding list of detected harmonics is presented 
in Table O It tabulates the multimode solution to the 
RV data set from SigSpec. Analytical 1-a uncertainties 
are ev aluated according to iMontgomerv fe Odonoghuel 
(1999); those in Table [U are 4-cr uncertainties. After 
prewhitening, the rms of the residual is 0.94 km s~ : and 
csig = 15.18. This indicates that the probability that 
the presented frequency solution could be generated by 
noise is 1 in 10 cslg . Figure [2}d plots the residuals after 
prewhitening. As shown, the residuals are small but not 
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TABLE 2 

RlGEL PERIODOCITIES FOUND FROM THE RADIAL VELOCITY STUDY. 



ID 




Frequency 




Amplitude 


(snr) rv 


sig 


r.m.s. 






[d- 1 ] 






[km s 


J 






[km s _1 ] 


fl 





01523 


± U 


,00002 


1 


128 


± 


0.027 


15.9 


105.9 


2.189 
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perfectly featureless. To avoid misidentification of true 
frequencies from aliases, we conservatively discontinued 
the prcwhitening at this step. 

4.3. Comparing Results of SigSpec and Peri odO^ 

A straightforward prewhitening of the data with the 
same weights used as above with Period04 starts with 
the same results as in SigSpec for the first few har- 
monics, but then shows some differences. Instead, we 
imported the frequencies from the SigSpec list, and 
prewhitened sequentially. The maximum difference of 
0.01 km s" 1 between the calculated amplitudes occurs 
in /19. The final rms of the residual is 0.92 km s" 1 in 
Perod04 and 0.94 km s _1 in SigSpec. The signal-to- 
noise-ratio within a box of 1.0 d _1 around fig is 4.70. 

4.4. Flux Amplitudes 

Independent prewhitening of the MOST dataset re- 
sults in tens of modes that have no corresponding coun- 
terparts in the frequency list of Table [5J This is quite 
expected since the baseline of MOST observations (~28 
days) is shorter than some of the periods found for the 
star indicated from the analysis of the RV data. To re- 
solve this, we used the "fixed" frequencies derived from 
the RV analysis (Tablc[2]) and employed Period04 to de- 
termine the possible corresponding light amplitudes in 
the MOST data. We searched for stable amplitudes and 
phases (an additional harmonic at the orbital frequency 
of the satellite / or b = 14.19827 d _1 was also subtracted). 
Unfortunately, due to the presence of long-period modes 
comparable to (and longer than) the MOST photometry 
baseline, no converging solution could be achieved. 

4.5. Frequency Distribution of Multimode Solution 

According to the asymptot ic theory of non-radial stel- 
lar oscillations (Tassoul 1980), high order p- and g- modes 
presen t regularity in fre quencies and periods, respec- 
tively (|Aerts et al.ll2010"al) . In the case of Rigel and sim- 
ilar classes of stars, which are believed to be unstable 
against g-modes, this regularity in period can bring a 



wealth of information about how s tellar material is mix- 
ing in the stars' deep interior (see iDegroote et aLll2010l 
for an example). However, as a prerequisite for applying 
this technique, identification of polar and azimuthal de- 
grees (£, m) for each individual mode is necessary. Unfor- 
tunately, the number of secured frequencies is not large 
enough for a statistical study of any possible period spac- 
ings. Moreover, the mode identification is beyond the 
scope of this paper. 

However, after extraction of the frequencies with good 
precision, it is worthwhile to demonstrate how the modes 
are distributed according to their corresponding frequen- 
cies and amplitudes. For this purpose, the left panel in 
Figure [5] shows the distribution of detected frequencies 
between [0, 1] d _1 and their corresponding RV ampli- 
tudes. Except fie and /19, the rest of the modes have 
frequencies below 0.4 d _1 . The right panel of Figure [3] 
shows the corresponding period distribution. 

It is possible that the rotation period of Rigel could 
be detected i n the frequency analyses of the RV data. 
IPapics et al.l (|2011[ ) found the rotational modulation as 
the only explanation for the variability in the B0. 5 
IV C0R0T double-line spectroscopic binary HD 51756. 
Rigel has t wo modern, self-consistent spectroscopic v sin i 
measures (jPrzvbilla et al.ll2010l ISimon-Di'az et al.ll2010D 
from which we adopt i>sini — 25 km s . This vshii 
is typical for a blue supergiant and indicates that an in- 
clination of i > 60° is likely (the inclination i is defined 
in the usual way as the angle of star's rotation axis rel- 
ative to our line of sight). Adopting both i = 60° and 
i = 90° with R = 78.9i? , yield P Iot (i = 60°) = 137 d 
(=0.0073 d- 1 ) and P mt (i = 90°) = 158 d (=0.0063 d" 1 ), 
respectively. The value of i > 60° is also in accord with 
the simple angular momentum conservation assumption 
from the expected w ro t = 240-300 km s _1 estimated for 
Rigel when it was a main-sequence O9/B0 V star. As 
shown in Table [2J there are no significant frequencies in 
the range of 0.006 - 0.007 d _1 . So, rotation effects can 
be neglected in our interpretation of frequencies. Yet, 
we are aware that multiple integers of the rotational fre- 
quency might show up in the frequency analysis, and 
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ulations induced by stellar spots are unlikely. 

4.6. Frequency Analysis of Equivalent Widths 

As panels c to e in Figure [2] clearly show, Si II A6347, 
Si II A6371, He I 6678 and Ha exhibit the large ampli- 
tude changes in EW. However, the binned EW measures 
have a very low duty cycle and DFT analyses do not 
lead to any significant frequencies above the 4-cr noise 
threshold. There is only one mode detected for the Si II 
A6347 line with frequency, amplitude and SNR of 0.0336 
d _1 , 14.71 mA and 4.24, respectively. Similarly for Si 
II A6371 line, the only dominant peak in the Fourier 
power spectrum corresponds to the frequency, amplitude 
and SNR of 0.0336 d" 1 , 11.32 mA and 3.46, respectively. 
Therefore, not only do the frequencies from EWs show 
no apparent correspondences with any of the frequency 
entries in Tabled the first two smallest frequencies, i.e. 
fs and /i, do not appear in the DFT power spectrum 
of EWs. Hence, with the current time series of spectra 
at our disposal, we cannot associate the low frequency 
range of RV variations with EW variations. 

5. SUMMARY 

This is the first time that simultaneous space-based 
photometry and time-resolved medium-resolution optical 
spectroscopy of a blue supcrgiant star are presented. Al- 
though the short baseline of MOST photometry did not 
enable us to derive flux amplitudes, the long baseline of 
RV monitoring and the DFT analysis revealed 19 signif- 
icant modes above SNR=4.6. These have been shown to 
arise from non-radial puls ations by | Kaufer et aT] (|1997f ) 
with gravity-mode nature (jLefever et al.ll2007f ). 

The question of the degree of regularity in light and RV 
variability in BA supergiants remains unanswered until 
an uninterrupted long time-baseline space-based observa- 
tion campaign of these objects is carried out. It is pos- 
sible that a degree of semi-regularity might exist in the 
current RV data set. But, there is no definite evidence 
for this, and thus it was not included in our frequency 
extraction procedure. Therefore, we caution that the pc- 
riodocitics found in our analysis might be affected by this 



ried out over several months by the BRITE-Constellation 
Mission dKuschnig et al.|[2009D . Observations like these 
should result in better defined frequencies. 

The temporal modulations in the equivalent widths of 
the two Si II lines show a completely different pattern 
from those of He I A6678 and Ha line, which indicate that 
equivalent widths formed at different optical depths t\ 
are uncorrelated. This demonstrates that there exists a 
pulsationally induced non-radial velocity gradient across 
different optical depths from which the lin es are formed 
(|Aerts et all 120091: ISimon-Di'az et al.ll2010l ). As a result, 
the response of the dilute extended atmosphere of Rigcl 
(which is reminiscent of other BA supergiants) to the 
pulsational waves emerging from the interior of the star 
is very complex, and has a high potential to probe a 
wide range of depths where the spectral lines are formed. 
As shown in Figure HJi, the equivalent width of the Ha 
feature shows the largest systematic variations of any of 
the spectral features measured. Because the Ha feature 
is also affected by mass outflows/inflows, and winds, it is 
a difficult line to interpret. It is, therefore, not surprising 
that its behavior is different than the other spectral lines 
studied. A full asteroseismic study of the Rigcl is planned 
in a subsequent paper. 
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